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It has long been recognized that observation of the 7 rays originating from nuclear deexcitation can 
be exploited to identify neutral-current neutrino-nucleus interactions in water-Cherenkov detectors. 
We report the results of a calculation of the neutrino- and antineutrino-induced 7-ray production 
cross section for the oxygen target. Our analysis is focused on the kinematical region of neutrino 
energy larger than ~200 MeV, in which a single-nucleon knockout is known to be the dominant 
reaction mechanism. The numerical results have been obtained using for the first time a realistic 
model of the target spectral function, extensively tested against electron-nucleus scattering data. 
We find that at a neutrino energy of 600 MeV the fraction of neutral-current interactions leading 
to emission of 7 rays of energy larger than 6 MeV is ~41%, and that the contribution of the p 3 / 2 
state is overwhelming. 

PACS numbers: 25.30.Pt, 21.10.Pc, 23.20.Lv 



The observation of 7 rays originating from nuclear de- 
excitation can be exploited to identify neutral-current 
(NC) neutrino-nucleus interactions in a broad energy 
range. The authors of Ref. [l| first suggested to use 
this signal to detect supernova neutrinos, the average en- 
ergy of which is ~25 MeV. Interactions of atmospheric 
neutrinos, with energies extending to the GeV region, 
can also lead to transitions to excited nuclear states de- 
caying through 7-ray emission, possibly associated with 
a hadronic cascade [2|. 

Neutrons, while providing ~50% of NC events, do not 
emit Chercnkov light. As a consequence, the availabil- 
ity of an alternative signal allowing one to identify NC 
interactions is very important. Events with 7 rays of 
energy above the observational threshold of 5 MeV can 
be detected in a water-Cherenkov detector, like Super- 
Kamiokande, and contribute up to ~5% of the total event 
number 0, independent of neutrino oscillations. Note 
that in water ~90% (16 out of 18) of the NC interactions 
take place in oxygen. 

Following the pioneering studies of nuclear excitations 
by neutral weak currents of Refs. p], @, theoretical cal- 
culations of the cross section of 7-ray production from 
NC neutrino-oxygen interactions have been carried out 
in the neutrino energy range E v ~ 10-500 MeV [j], 0, [![■ 
These studies took into account 7 rays originating from 
the inelastic processes v + ^O v 1 + l %0*, in which 
the oxygen nucleus is mainly excited to resonances ly- 
ing above particle emission threshold. These states then 
decay to either p + ^N* or n + 1 |0 , and the residual 
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nuclei, left in excited particle-bound states, decay in turn 
emitting 7 rays in the 5-10 MeV region. 

At low energy, elastic scattering and inelastic excita- 
tion of discrete nuclear states provide the main contribu- 
tion to the neutrino-nucleus cross section. However, at 
E v > 200 MeV the cross section associated with these 
processes tends to saturate, and quasielastic (QE) nu- 
cleon knockout becomes the dominant reaction mecha- 
nism. If the residual nucleus is left in an excited state, 
these processes can also lead to 7-ray emission. The K2K 
Collaboration reported the observation of 7 rays from 
nuclear deexcitation following NC neutrino-oxygen inter- 
actions at E v ~ 1.3 GeV in the 1-kton water-Cherenkov 
detector 0] . The number of events and the visible energy 
are qualitatively consistent with those expected from 6- 
MeV 7-ray production in NC QE neutrino-oxygen inter- 
actions. 

In the QE regime, neutrino-nucleus scattering reduces 
to the incoherent sum of elementary scattering processes 
involving individual nuclcons, the energy and momentum 
of which are distributed according to the target spec- 
tral function [l(j ■ A schematic representation of NC QE 
neutrino-nucleus scattering is given in Fig. [TJ where the 
dashed line represents the threshold for nucleon emission 
in the continuum. 

In this Letter, we discuss the emission of 7 rays arising 
from the decay of the residual nuclei of the reactions v + 
^O -> p+p + ^N* or u + ^O -> u + n+^O*, the cross 
sections of which have been computed using a realistic 
model of the oxygen spectral function. 

Note that, due to the strong energy- momentum corre- 
lation exhibited by the nuclear spectral function, large 
excitation energies of the residual system are associ- 
ated with large momenta of the knocked out nucleon. 
As nucleons occupying shell-model states have a van- 
ishingly small probability of carrying momentum larger 
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FIG. 1. (color online). Schematic representation of neutral- 
current neutrino scattering off oxygen. 
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FIG. 2. Low-lying excited levels of the residual nuclei pro- 
duced in 1 |0(^, v' N) scattering. Energies are measured with 
respect to the ground state. 



than ^250 MeV [111, knockout of these nucleons predom- 
inantly leaves the residual system in a bound state. 

In our approach, the cross section of 7-ray production 
following a NC QE interaction, a 1 , is written in the form 

<t 7 = cr(v + 1 |0 -^-v + j + Y + N) (1) 

= E a ( v + W s° v + X » + N ) Br ( X « 7 + Y), 

a 

where iV is the knocked out nucleon, X a denotes the 
residual nucleus in the state a, and Y is the system re- 
sulting from the electromagnetic decay of X a , e.g. 



15 7 N, 



15q 



or gC +p [12Hl4j . The energy spectrum 
of the states of the residual nuclei is schematically illus- 
trated in Fig. H 

According to the shell model, nuclear dynamics can 
be described by a mean field. In the simplest imple- 
mentation of this model, protons in the 1 gO nucleus oc- 
cupy three states, lpi/2, I.P3/2! an d lsi/2, with removal 
energy 12.1, 18.4, and ~42 MeV, respectively [l5|-[l7|. 
The neutron levels exhibit the same pattern, see Fig. [TJ 
but are more deeply bound by 3.54 MeV [14]]. Since be- 
low nucleon-emission threshold the deexcitation process 
is governed only by energy differences, the proton and 
neutron holes yield photons of very similar energy, the 



differences being as small as ~0.1 MeV (see Fig. [2]). 

The calculation of the NC QE cross section, a(v + 
1 |0 — > v + X a + N), has been performed within the ap- 
proach discussed in Refs. for the case of charged- 
current (CC) interactions, whereas the branc hing ratios 
Br(X a 7 + Y) have been taken from Refs. (TJ HI- 

Following Refs. [H, El, we write the NC QE cross 
section in the form 
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where E N = ^AP + p 2 , M being the nucleon mass, 
da U N I d£ldE' v denotes the elementary neutrino-nucleon 
cross section and the spectral function Pjy(p,E) yields 
the probability of removing a nucleon of momentum p 
from the target leaving the residual nucleus with energy 
E + E — M, Eq being the target ground-state energy. 

In the nuclear shell model, nucleons occupy single- 
particle states <p a with binding energy — E a {E a > 0). 
As a consequence, knockout of a target nucleon leaves 
the residual system in a bound state, and the spectral 
function can be conveniently written in the form 



Pn(p,E) 



E 

a 6{F} 



n a \Mp)\ 2 UE-E a ), (3) 



where 4> a (p) is the momentum-space wave function as- 
sociated with the ath shell model state and the sum is 
extended to all occupied states belonging to the Fermi sea 
{F}. The occupation probability n a < 1 and the (unit- 
normalized) function f a (E — E a ) , describing the energy 
width of the ath state, account for the effects of nucleon- 
nucleon (NN) correlations, not included in the mean-field 
picture. In the absence of correlations, n a — > 1 and 
f a {E - E a ) -> 5{E - E a ). 

Precise measurements of the coincidence (e, e'p) cross 
section, yielding direct access to the target spectral func- 
tion, have provided unambiguous evidence of deviations 
from the mean-field scenario, leading to significant de- 
pletion of the single- particle states fllMTl ] . The data at 
large missing momentum and large missing energy [i.e. 
large |p| and large E in Eq. collected at Jefferson 
Lab by the JLAB E97-006 Collaboration, indicate that 
NN correlations push ^20% of the total strength to con- 
tinuum states outside the Fermi sea (2lj . 

A realistic model of the proton spectral function of 
oxygen has been obtained within the local density ap- 
proximation (LDA) , combining the experimental data of 
Ref. 15( with the results of theoretical calculations of 
the correlation contribution in uniform nuclear matter 
at different densities [l8|, [22J. The results reported in 
Ref. [3 show that the LDA spectral function provides 
an accurate description of the inclusive electron-oxygen 
cross sections at beam energies around 1 GeV. In ad- 
dition, it predicts a nucleon momentum distribution in 



agreement with that obtained from the data of Ref. [21 

As pointed out by the authors of Ref. [23}, nucleon- 
knockout experiments measure spectroscopic strengths, 
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FIG. 3. (color online). Distribution of removal energy of 
protons and neutrons in x |0. 



TABLE I. Spectroscopic strengths of the 1 |0 hole states and 
their branching ratios for deexcitation by the _E 7 > 6 MeV 
photon emission. 



a 


Pl/2 


P3/2 


Sl/2 




0.632 


0.703 


0.422 


Br(X a 7 + Y) 


0% 


100% 


16 ± 1% 



not occupation probabilities. Spectroscopic strengths are 
given by the area below the sharp peaks observed in 
the missing-energy spectra, corresponding to knockout 
of a nucleoli occupying one of the shell-model states, cor- 
rected to take into account final-state interactions. On 
the other hand, occupation probabilities include contri- 
butions corresponding to larger removal energy, arising 
from mixing of the one-hole state with more complex fi- 
nal states [23f . 

Figure [3] shows the energy distribution obtained 
from momentum integration of the spectral function of 
Refs. [IHlIIj]. It clearly appears that, unlike the pi/ 2 and 
P3/2 states, the S1/2 state is spread out over a broad en- 
ergy range, and can hardly be treated as a single-particle 
state. 

The P1/2, P3/2 and S1/2 spectroscopic strengths have 
been computed by integrating the oxygen spectral func- 
tion of Refs. 0, [13 over the energy ranges 11.0 < E < 
14.0 MeV, 17.25 < E < 22.75 MeV, and 22.75 < E < 
62.25 MeV, respectively. Dividing these numbers by the 
degeneracy of the shell-model states, one obtains the 
quantities S a listed in Table HI The same spectroscopic 
strengths have been used for protons and neutrons. 

Our results turn out to be very close to those ex- 
tracted from the high resolution measurement carried 
out at NIKHEF-K [lg. For example, the pi/ 2 (P3/2) 
strength collected in the same energy range is reported 
to be 0.630±0.034 (0.676±0.037). 

The uncertainty in the determination of S a is mainly 
due to the choice of the shell- model wave functions and to 
the treatment of final-state interactions of the knocked- 
out proton. The authors of Ref. [l|| quote an overall 
systematic uncertainty of 5.4%. 



The elementary neutrino-nucleon cross section of 
Eq. ^ can be written in the form 
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where E' N = y M 2 + p' 2 . The leptonic and hadronic 
v , are given by 

(5) 



tensor, and W 1 *" , are given by 

L^ LV = 2(fc^fc„ + k u k^ — g^yk ■ k — %£^ va sik k ^) 
and 

+ q^q v Wi + (p^q v + p"r)W 5 , 



(6) 



with p — (En, p) and q = (£»,k — k'). As in the case of 
CC QE scattering [3, the structure functions Wi can be 
written in terms of the nuclcon form factors according to 
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with r = — q 2 /(AM 2 ). Note that, in the above equa- 
tions, the electromagnetic and charged-current nucleon 
form factors {Ff (i — 1,2), Fa, F p } arc replaced by the 
ones appropriate to describe NC interactions @, [2J, [25[ 
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where the upper (lower) sign corresponds to proton (neu- 
tron) form factors, 0w is the weak mixing angle, is 
the pion mass, g A = —1.2673, and the strange quark con- 
tribution is set to As = -0.08 The form factors Ff 
and F 2 can be expressed in terms of the measured Sachs 
form factors G^ and as 
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(9) 



In this Letter, we use the state-of-the-art parametrization 



of Gg and G% of Ref. [27 



The branching ratios Bi~(X a —> 7 + Y), necessary to 
calculate the cross section tr 7 according to Eq. (JT]), arc 
collected in Table U [lH |2(j. In the case of the Pi/ 2 - 
proton (neutron) knockout, the residual nucleus is 
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FIG. 4. (color online). Cross section for 7-ray production fol- 
lowing NC QE interaction of neutrino (solid line) and antineu- 
trino (long-dashed line) compared to the NC QE cross section 
of neutrino (dotted line) and antineutrino (short-dashed line). 
Only the photons of energy larger that 6 MeV are considered. 



( |0) produced in its ground state. Hence, no 7 rays are 
produced. As the p 3 / 2 -proton (neutron) hole lies below 
the nucleon-emission threshold, 10.21 MeV (7.30 MeV), 
it always deexcites through photon emission with half-life 
0.146 ± 0.008 fs (less than 1.74 fs) QJ]. When a proton 
(neutron) is knocked out from the deepest S1/2 shell, the 
excitation energy is high enough for many deexcitation 
channels to open, of which only two, ^C+p and ^N + n 
( 1 |C + p an d 1 p C + a), yield photons of energy higher 



than 6 MeV 12J (see Fig. (2). The theoretical estimate of 
the branching ratio for these processes [HI, [2(| , being in 
total 16%, turns out to be in good agreement with the 
value 15.6 ± 1.3±5;^% extracted from the ^Oip^p) 1 ^ 
experiment El 48 carried out at the Research Center for 
Nuclear Physics (RCNP) of the Osaka University [28] . 

Although the Sj/2 contribution to cr 7 is an increasing 
function of neutrino energy, it saturates at ^0.4 GeV and 
remains less than 4.5% for E v < 5 GeV. Compared to 
the contribution, it is suppressed by the low branch- 
ing ratio for deexcitation through 7 emission, the lower 
strength and degeneracy, as well as the larger removal 
energy, that makes the NC QE cross section smaller. 

Note that the formalism presented in this Letter also 
applies to antincutrino-induccd 7-ray production, the 
only difference being the sign of the last term in Eq. . 

In Fig.[4]our results for the neutrino- and antincutrino- 
induced 7-ray production cross section are compared to 
the neutrino and antineutrino NC QE cross sections. The 
error bands show the uncertainties arising form the deter- 
mination of the spectroscopic strengths (5.4%), the treat- 
ment of Pauli blocking (1%), and the branching ratio of 
the S1/2 state (1%). The <7 7 's dependence on neutrino 
energy is very similar, although not identical, to that 
of the NC QE cross section. The discrepancy arises from 



difference between the average removal energy associated 
with the whole spectral function and the energy of the 
p 3 / 2 shell, yielding the overwhelming contribution to a 1 . 
The neutrino-induced 7-production cross section reaches 
its maximum at E v ~ 1.9 GeV and is slowly decreasing 
at larger energies. On the other hand, the corresponding 
antineutrino cross section is an increasing function of E v . 

As the axial mass enters a 1 only through the elemen- 
tary cross section the ratio ct 7 /<7ncj ^nc being the 
NC QE cross section, is largely independent of Ma- For 
example, applying Ma = 1-39 GeV (1.03 GeV) instead of 
the value 1.2 GeV used in this Letter [2^|, changes the ra- 
tio cr 7 /cr N c by less than 0.4% (0.3%) when E v < 5 GeV. 

The mechanism of 7-ray production through nuclear 
deexcitation is the same for NC and CC processes. 
Therefore, the fraction of neutrino interactions followed 
by 7-ray emission in the two cases is determined by the 
same factors. In the case of CC QE scattering, the max- 
imum value of the energy transfer is reduced by the non- 
vanishing mass of the charged lepton, and therefore the 
contribution of the p 3 / 2 knockout to the total cross sec- 
tion is somewhat more significant. However, this effect 
is already small at E v = 475 MeV, as the fraction of CC 
interactions emitting photons is higher than that of NC 
interactions by only 1%, and becomes even smaller with 
increasing neutrino energy. 

In conclusion, we have computed the neutrino and 
antineutrino neutral-current cross sections, focusing on 
the kinematical region in which single-nucleon knockout 
dominates. In this region the average of neutrino and an- 
tineutrino cross sections obtained from our approach is 
much larger than the corresponding result of Ref . Q . For 
example, at E v = 0.5 (1.0) GeV our average cross section 
exceeds the one reported in Ref. [1] by a factor ~10 (15). 

The NC cross sections have been used to compute the 
7-ray production cross sections. Considering photons of 
energy larger than 6 MeV, we find that the p 3 / 2 state 
provides the overwhelming contribution, and that the ra- 
tio <7 7 /ctnCj exhibiting a significant energy-dependence at 
E v < 1 GeV, is ~41% at E v = 600 MeV. 

Our results, obtained using a realistic model of the 
target spectral function, provide an accurate estimate of 
a signal that can be exploited to identify neutral-current 
events in water-Cherenkov detectors. 
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